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Introduction 

In  many  aerodynamic  experiments,  the  mathematical  model  of  measurements  is  the 
following  equation 

fn  ~  ^  ^ n  Vn  5  n  —  2,  •••  5  ( 1 ) 

where  6  is  the  unknown  vector  of  dimension  M,  kn  is  the  “scattering”  function  of  the 

measuring  system  (string-vector  of  dimension  M),  r]n  is  the  noise  of  measurements  with  zero 
average  value  and  dispersion  a 2 .  The  problem  of  interpretation  of  such  measurements 
consists  in  construction  of  an  estimation  0)  n>  for  the  vector  0  on  the  basis  of  measurements 
fi,i  =  1,  2,  ...,  n  .  This  problem  is  solved  using  recurrent  algorithms,  in  which  an  update  of  the 

estimation  0in^ 1  (constructed  on  the  previous  step)  is  implemented  at  each  new  measurement 

/,, .  The  recurrent  algorithm  is  described  by  equations 
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with  a  start  point  A  zero  vector  is  often  assumed  as  a  vector  6^°' .  A  matrix  of 

dimension  M  xM  is  set  by  the  expression  a-1  •  /  ,  where  a  is  the  regularization  parameter 

and  /  is  the  unit  matrix.  It  is  evident  that  the  error  of  the  estimate  6^”*  depends  on  the  step 
number  n  .  However,  the  choice  of  the  recurrent  algorithm  stop  moment  is  considered  in  the 
literature  rather  insufficiently.  Therefore  a  question  arises  “How  to  choose  the  number  n  for 
which  the  norm  of  the  solution  error  was  minimum  or  recurrent  algorithm  meets  certain 
requirements”.  To  answer  this  question,  the  accuracy  characteristics  of  the  recurrent  algorithm 
(2),  (3)  used  to  solve  the  problem  of  algorithm  synthesis  are  introduced.  It  should  be  noted, 
that  the  solution  of  this  problem  allows  one  to  find  the  least  number  of  the  measurements, 

which  is  necessary  to  obtain  the  estimate  #*-"  *  with  the  required  accuracy. 

Accuracy  characteristics  of  the  recurrent  algorithm 
At  each  step,  the  vector  of  the  estimate  error  is  defined  by  the  expression 
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where  8  is  the  “exact”  vector  of  unknown  parameters.  The  vector  i'*"1  may  be  written  as  a 
sum  of  two  components  s ^  =  ^n]  +  b l'"'> ,  where  £*"■  =  8U !)  -  8 ("-  is  the  vector  of  the  random 

error  caused  by  “transferring”  to  the  solution  6[n)  the  random  error  of  initial  data; 
bin>  =  8{" ■  -6  is  the  vector  of  the  systematic  error.  The  vector  8 (,4  is  the  solution  of  system 
(2),  (3)  for  exact  values  /„  =  k„  -  8  .  Let  us  consider  the  calculation  of  these  vectors. 

We  will  enter  a  rectangular  matrix  B(n)  of  dimension  M  xM  involved  in  the  expression 

b(n)  =  B(n)  ■  (8(0)  -8)  =  B(n)b(0) , 


where  b(0)  =6»(0)  -8  is  the  vector  of  “initial”  bias.  We  will  call  B(n)  a  bias  matrix. 

Statement  1.  The  matrix  B(n)  for  the  recurrent  algorithm  (2),  (3)  is  defined  by  the 
expression 


n— 1 

B(n)  =  11 

7=0 

and  allows  for  recurrent  calculation: 
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■  B(n  -1) ,  n  =  1,2,... 


where  .5(0)  =  / . 

We  will  take  the  correlation  matrix  Vr(n)  =  M 


as  a  characteristic  for  the 


vector  of  a  random  error  £pn’ ,  where  M[  ]  is  the  operator  of  mathematical  expectation 
with  respect  to  the  noise  distribution. 

Statement  2.  The  correlation  matrix  Vr(n)  allows  the  recurrent  calculation: 


Vf(n)  =  (I -G  (n)  ■kn)-V^{n-\)-(I-G  (n)  ■  kn  y  +  of,  ■  G  (n)  ■  Gr  («), 


(4) 


Vg(0)  =  0,  «  =  1,  2,..., 

where 

P{n~X)kl 

G  («)  = - — - . 

knP^n  ^ kn  +  <T2 

The  introduced  bias  matrix  B(n)  and  correlation  matrix  f>(«)  characterize  the 

systematic  and  random  errors  sufficiently.  However,  in  order  to  simplify  the  synthesis,  we 
should  introduce  some  additional  characteristics  of  these  matrices. 

It  can  be  easily  shown  that  the  root-mean-square  error  (rms)  of  the  estimate  8 *-")  is 
defined  by  the  formula 

A 2(/i)  =  M  [||  8{n)-8  ||2  ]  =  Z,(,,)V,!)  +Sp[V^{n)],  (5) 

where  Sp  [  ]  is  the  matrix  trace.  Let  us  consider  each  addend  in  (5)  separately. 
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In  practice  it  is  impossible  to  calculate  the  vector  M") ,  because  we  do  not  know  the 
“initial”  bias,  that  is,  the  vector  b^f  Therefore  we  shall  define  the  scalar  characteristic 
Ub(n)  [1],  defined  as 
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where  \mx\  is  the  unit  vector,  bl^  is  the  bias  vector 
bias,  that  is,  b t0)  =  \mx\  •  We  can  show  that 


1 


M  f  M 


=  £{*(«)} 
1=1  v  ./=| 


stipulated  by  a  unit  vector  of  initial 
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The  following  limits  are  valid: 


lim  Ub(n)  =  1,  lim  Ub{n )  =  0  . 

n  — >0  n  — >oo 

Let  us  consider  the  second  addend  in  (5).  We  will  determine  the  scalar  characteristic 
Ug(ri)  [1],  defined  as 

that  can  be  interpreted  as  a  coefficient  of  transference  of  noise  variance  to  the  quantity 

M[||  dn)  |2]  =  S^(n)]. 

As  it  follows  from  (4),  the  algorithm  for  calculating  Ug(n)  can  be  defined: 


U£ (n)  =  \sp[( I  -  Ga (n)  -kn  )  ■  V^(n  - 1)  •  ( I  -  Ga{n)  ■  kn  f]  +  Sp [ Ga ( n )  •  GTa(n) ]  . 

Using  the  characteristics  Ub{n),  U b (n) ,  we  can  approximate  the  MSE  of  estimation 
by  the  expression 

ii  m2 

A2  («)  a  Ub(n )  •  6(0)  p  +  (n)  ■  cr2  .  (6) 

The  values  U b  ( n )  and  U ^  ( n )  determine  systematic  and  random  errors  too  frill  and  they 
maybe  called  the  accuracy  characteristics  of  the  recurrent  algorithm  (2),  (3). 

Let  us  consider  the  results  of  one  computing  experiment.  The  vector  6 ,  having  30 
projections,  is  measured  using  (1).  The  number  of  measurements  N  is  equal  to  80.  The 
condition  number  of  a  the  measurement  matrix  is  equal  2-107.  The  relative  level  of 
measurement  noise  is  equal  0.05.  Figures  1  and  2,  respectively,  demonstrate  Ub,U%  as 
functions  of  the  iteration  number  n  .  They  show  the  plots  for  two  regularization  parameters: 
a  =  10~*  (plot  1)  and  a  =10  7  (plot  2).  The  behaviour  of  Ub{n)  and  U < ( n )  demonstrate  a 
known  inconsistency  between  the  stability  and  the  resolution  of  the  algorithm.  With  increase 
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Fig.  1  Fig.  2 

of  the  iteration  number  n  the  resolution  is  improved ,  but  the  stability  of  the  estimation  0(n'1  to 
measurement  noise  is  worsened. 

It  is  evident,  that  calculating  these  characteristics  we  can  perform  an  accuracy  analysis  of 
the  recurrent  algorithm  (2),  (3).  For  this  purpose  in  formula  (6),  instead  of  an  unknown 

2 

quantity  6<0)  ,  its  estimation  b^ax  is  substituted.  To  illustrate  this  possibility,  the  plots  of 

the  functions 

At(n)=Ub(n)-bZmK+Uf(n)-cjZ;  (7) 


8{n) 
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calculated  for  the  previously  described  measurements  (see  Fig.  I  and  Fig. 2)  are  presented  in  the 
Fig.  3.  The  value  of  6,„ax  was  30  %  more  than  the  value  h<0)  .  From  Fig.  3  it  is  visible  that 

the  minimum  point  of  the  function  A;(«)  (random  error  of  the  estimation  (fn') )  is  in  a 
neighbourhood  of  the  minimum  point  of  the  function  Af(n)  (expectation  of  the  random  error 
of  estimation),  and  therefore,  A ^(n)  can  be  used  to  analyze  the  recurrent  algorithm  (2),  (3). 

An  important  property  of  A \{n) ,  Ub(n),Ug(n)  is  the  possibility  of  their  “a  priori” 

calculation,  since  the  measurements  f  are  not  involved  in  calculation  of  the  accuracy 
characteristics.  This  allows  us  to  determine  the  iteration  number  proceeding  from  the  required 
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Fig.  3 

values  of  these  characteristics,  i.e.  perform  synthesis  of  the  recurrent  algorithm.  For  this 
purpose,  we  will  formulate  the  following  variational  problems : 

Problem  A 

min  A ](ri)- 

Problem  B 

mmUb{n)  under  constraint  Up(n)  <  C/?ax  . 

«>0  5 

Problem  C 

mmUp{n)  under  constraint  Ub^UV1™. 

n>  0  b 

Let  us  note  the  following: 

-  The  solution  of  problem  A  (denote  it  as  )  is  an  estimate  of  the  iteration  number  at 
which  the  rms  error  of  the  estimate  0^  is  minimum,  but  as  a  priori  information  we  must 
specify  the  values  of  ||z/0'  ||  and  .  As  a  rule,  these  values  are  unknown.  Therefore,  we  shall 

construct  a  recurrent  algorithm,  which  will  be  optimum  both  on  classes  of  vectors  6  and 
measurement  noise  satisfying  the  following  limitations: 

I  b{0)  f  <  b ;  al  <  i  .  (8) 

Substituting  Z>^ax ,  cr^ax  *nt0  (7)  and  solving  the  problem  A  we  find  the  n ^ ,  that  is,  the 

number  of  steps  minimizing  the  rms  error  of  the  estimate  d^1'  on  classes  (8). 

-  The  solution  of  the  problem  B  (denote  it  by  ng )  minimizes  the  systematic  error  with 
the  required  stability  to  the  measurement  noise  . 

—  The  solution  of  the  problem  C  (denote  it  by  np )  minimizes  the  random  error  with 
the  required  algorithm  resolution. 
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-  Taking  into  account  the  monotone  behavior  of  the  characteristics  U b(ri),U p{ri)  ,  the 
solution  of  the  variational  problem  B  is  determined  as  the  root  of  the  nonlinear  equation 

U4{nB)  =  Uf™, 

and  the  solution  of  the  variational  problem  C  -  as  the  root  of  the  nonlinear  equation 

Ub(nc)  =  C/,max  . 

The  sign  =  shows  that  the  solution  is  sought  among  the  integer  values  of  the  arguments 
Ub(n)  and  U^{n) . 


Conclusion 

The  proposed  accuracy  characteristics  of  the  recurrent  algorithm  (2), (3)  allow  us  to 
perform  not  only  its  analysis  but  also  the  synthesis,  i.e.,  to  choose  the  iteration  number  of  the 
recurrent  algorithm  stopping,  at  which  the  desired  accuracy  characteristics  are  ensured. 
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